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Abstract 

We study two-dimensional chaotic standard maps coupled along the 
edges of scale-free trees and tree-like subgraph (4-star) with a non- 
symplectic coupling and time delay between the nodes. Apart from 
the chaotic and regular 2-periodic motion, the coupled map system ex- 
hibits variety of dynamical effects in a wide range of coupling strengths. 
This includes dynamical localization, emergent periodicity, and ap- 
pearance of strange non-chaotic attractors. Near the strange attrac- 
tors we find long-range correlations in the intervals of return-times to 
specified parts of the phase space. We substantiate the analysis with 
the finite-time Lyapunov stability. We also give some quantitative ev- 
idence of how the small-scale dynamics at 4-star motifs participates 
in the genesis of the collective behavior at the whole network. 

Keywords: Self-organized criticality (Theory); Network dynamics; New ap- 
plications of statistical mechanics; 

1 Introduction 

Network Dynamics, Diffusion & Chaos. Complex dynamical systems can 
be efficiently modeled by network topologies where the edges represent the 
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interactions between the dynamical variables attached to the network nodes 
[H [21 [3]. Various dynamical processes on the networks have been studied 
in which the temporal fluctuations in the dynamical variables constrained 
by the network structure are subject to the inter- node interactions and/or 
driving by the external noise. The interplay between the network topology 
and emergent dynamical behavior can be readily demonstrated in diffusive 
processes on scale-free networks [3J. Of particular importance for the dy- 
namics are the networks with well developed modular structures, topological 
subgraphs, which determine their own time scales [I]. Such structures are of- 
ten present in biological systems, in particular, in gene interaction networks, 
where structural modules, or dynamical motifs [5], represent gene functional 
units mutually connected through the network. Although the cooperativ- 
ity in gene functions can be easily demonstrated, e.g., with the appropriate 
statistical analysis [6] of the empirical gene-expression data, both the dynam- 
ics of gene motifs and the structure of gene networks remains a challenging 
problem. Recently, an application of discrete-time dynamics has been sug- 
gested for modeling gene interactions [H El E] • This approach brings a new 
prospective, in particular with respect to the advanced stability analysis and 
an intriguing possibility of a chaotic behavior. 

The coupled chaotic maps on regular networks [TUl EH CE2], structured 
small graphs [131 HH EH] an d scale-free networks [161 EH EH ESI [20] are cur- 
rently the subject of an intensive study. Many of the known one- dimensional 
chaotic maps were studied with the emphasis on the effects of coupling 
strengths and of the time delay [211 [221 [23] between the units. The central 
points in these studies are the phenomena related with the synchronization 
of spatially extended coupled maps [TH [151 ESI EHl l2H (221 [23] , nature of the 
transition to the synchronized state [TU] and coexistence of different types of 
attractors [13] . 

Generally, coupled chaotic maps are interesting dynamical systems in 
which the dynamic stability can be examined with respect to variations in: 

• Initial conditions (the trajectory divergence); 

• Parameters (coupling at a node, coupling between the units, geometry, 
time delay); 

• Drivings (external frequencies, endogenous driving). 

In this respect, two-dimensional coupled chaotic maps, which are much less 
studied so far, leave more room for the emergent dynamical effects [TTll2T)i rr4]. 
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In particular, in a system of standard maps with symplectic coupling on a 
full graphs trapping of the trajectories and anomalous diffusion was found 
[13] . In the coupling with time delay on scale- free tree the standard maps 
lead to characteristic dynamical patterns with robust statistical properties 
[20] . Type of the coupling between two degrees of freedom for a general one- 
dimensional map coupled along a linear chain was shown [IT] to affect the 
transition to the synchronized state. 

In this work we study the discrete versions of the two-dimensional peri- 
odically kicked oscillators, usually referred to as the Chirikov standard maps 
[23], coupled along their phase variables on large scale-free trees and their 
typical subgraph, the 4-star motif. As described in detail below, we impose 
a fixed time delay between coupled neighboring units (maps) thus rendering 
the coupled maps system non- symplectic. Our main purpose is to character- 
ize the dynamic behavior that emerges as a consequence of the inter-node 
coupling both at small scale structure - the 4-star motif, and at the large scale 
structure - the tree with N = 1000 nodes. Between the chaotic phase, which 
persists at weak coupling, and the fully synchronized phase with 2-periodic 
orbits, which dominates when nodes are strongly coupled, we find a whole 
range of couplings with a self-organized behaviors in which nodes influence 
each other. In this region the dynamical orbits are generally non-chaotic, 
with different patterns of collective dynamic behaviors, that we devote most 
of our attention to in this work. 

The paper is organized as follows. After defining our coupled map system 
in Section we report the qualitative and statistical features of the coupled 
maps in Section [3j In Section H] we investigate the dynamic stability of 
system's emergent orbits and describe behavior near the fractal attractors. 
Finally in Section we give a short summary and discussion of the results. 

2 The Graphs of Coupled Maps 

We consider the standard map [23] in its original unbounded version, with x 
being the phase (angle) and y being the action (momentum) variable defined 
on each node [i] of a network with iV nodes: 

/ x'[i] \ _ f x[i] + y[i] + e sm(27rx[i]) modi 
V S/'M / \ y[i] + £sm(27rx[i]) 
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This map has been extensively studied as an isolated system and its proper- 
ties are well known, in particular the ones regarding the phase space diffusion 
for e-values above the chaotic transition threshold (e c ~ 0.155). This is a 
symplectic map (a discrete version of a Hamiltonian system) which has a 
very intricate phase space structure with invariant KAM-torai persisting be- 
low the e c value [25] . 

We define our Coupled Map System (CMS) in accordance with the usually 
adopted scheme [17] as: (1 — fi) x (standard map-update) + fix (coupling). 
Also in analogy with the previous works on ID maps [IH [TTJ [26] , we consider 
the diffusive phase-coupling in the x- variable of Eq. (jTJ) . A one-step time- 
delay between the neighboring nodes is imposed, so that our CMS reads: 



Here, (') denotes the next iterate of Eq. (JTJ), n is the global discrete time, 
fi is the coupling strength, [i] indexes the nodes on the graph, ki is node 
degree and /Q denotes the neighborhood of the node [i]. The update of each 
node is the sum of a contribution given by its standard map-update (the ' 
part) plus a coupling contribution given by the sum of differences between the 
node's phase- value and the phase- values of neighboring nodes in the previous 
iteration, normalized by the node's degree. 

In contrast to the well-studied phase-coupled ID oscillators [TfJ [261 fl2j . 
our CMS defined by the Eqs. ([1][2]) is essentially a network of time-delayed 
phase-coupled ID oscillators with an additional coupling to the local mo- 
mentum variable at each node. From a continuous-time viewpoint, seeing 
x'(i,t) = x[i]'(t) as a continuous variable over a discrete space-time defined 
by the network, this coupling functional form can be seen as: 

*'(*, t + 1) = (1 - Ax®' + £ [£, 6X (./•/' - s[*T) + £ j6 K>b1 - *I>T) 
(i - im) x [%\' + f [v 2 x[^]' - z j& K t 9 MJy] - [(i - /*) + £(v 2 - hdt)] x'(i, t) 

(3) 

which represents a sort of diffusive phase-coupling on network. In compari- 
son to the diffusive CMS that are typically studied on chains (e.g. [TT]). our 
model involves 2D elements on a network with a time-delay, causing addi- 
tional dynamical effects. Note also that, in contrast to the study in Ref. [T4] . 
our coupling scheme defined by the Eq. (T5]) is non-symplectic and in addi- 
tion it includes a fixed time-delay. This type of coupling is motivated by a 
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phenomenological model of the node interaction coming from the oscillatory 
origin of the standard map. Another important feature of the CMS defined 
by Eq. ([2]) is its systematic inhibition of the diffusion along the action coordi- 
nate, which is made possible by the factor (1 — /i) in front of y[i]' (see section 
|3j), thus allowing more coherent interactions among the coupled nodes. 

Geometries. The CMS of Eq. (J2J) interact along the links of a network 
determined by its adjacency matrix C^. In this work we primarily focus 
on the coupling strength effects. Therefore, we consider tree-like structures 
in order to avoid possible loop effects on the time-delayed interactions. In 
particular, we study CMS defined by Eq. j2} on a large scale-free tree with 
N = 1000 nodes and on a typical 4-star motif, shown in Fig.[TJ Our previous 
study on trees [2U] suggests that the chaotic maps with the coupling in Eq. (j2J) 
stabilize through the localization of orbits into well defined bands in the 
action coordinate. These bands determine different clusters of periodically 
synchronized nodes, which make a statistically stable pattern on the tree 
with the characteristic distance d = 2 between two synchronized nodes [20] . 
Note that the 4-star is a typical motif which captures both the characteristic 
distance of synchronization and the tree-like nature of the topology. Here we 
systematically compare the emergent dynamic behavior of a 4-star and the 
whole tree in an extended range of couplings, where the CMS exhibits new 
dynamical phenomena, as we describe below. 

In Fig.d^i we show a typical situation of the tree at the threshold of 
localization as the coupling \i approaches a critical value jj, c ~ 0.021 from 
below. In colors are shown nodes which first achieve periodic orbits, in 
contrast to black nodes, at which the maps still remain chaotic. The figure 
shows that the localization of orbits on the tree starts from the nodes with 
least number of links. In Fig.[Tb we show the adjacency matrix and the 
situation with fully synchronized nodes that occurs at large coupling strength 
/i = 0.08. (For clarity a smaller scale-free network was shown.) An interesting 
feature is that the pattern of synchronized nodes shows the largest density 
at the top-right corner, representing nodes at the boundary of the graph, 
whereas the pattern of direct links between nodes has largest density at the 
origin (near the hub node). The results are collected from and averaged over 
many initial conditions. 

In the following two sections we describe different types of dynamical be- 
haviors obtained at different coupling strengths. Our particular emphasis is 
on the region 0.021 < \i < 0.08, i.e., between the above mentioned initial 
localization of orbits and the synchronization of nodes at strong couplings. 
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Figure 1: (a) Scale- free tree at the threshold of regular dynamics with a 
pattern of nodes (bright color) which first achieve regular orbits; Inset: 4- 
star motif used in this work, (b) Pattern of connections (yellow) and pattern 
of synchronized nodes (red) on a 100-node scale-free tree for large inter-node 
coupling defined by Eq. (j2|). 



We focus on the statistical properties of the collective motion and on the dy- 
namic stability of emergent orbits. Throughout this work we fix the coupling 
e = 0.9 between two components (x[i], y[i}) at each node and the time delay 
t = 1 between coupled nodes and we vary the inter-node coupling /i. The 
results are collected after an initial transient, typically 100000 steps for each 
trajectory, and averaged over many (typically 1000) initial conditions. In all 
runs we pick the initial condition for all the nodes of the network randomly 
with a uniform probability from the phase space subset (x, y) 6 [0, 1] x [— 1, 1]. 

3 Orbits Structure and Dynamical Regimes 

To study CMS on networks, the following types of emergent orbits (e.o.) can 
be considered: 

• Orbits of each individual node of the CMS, defined as 

(z[*k,y[*]n) re>no , » = 1, •• • , AT (4) 
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where no denotes the end of transients. For a selected node on the 
network we can in this way investigate in detail 2D discrete dynamics 
using well known study techniques of the dynamical systems. 

• Network-averaged emergent orbit (n.a.e.o.) defined as 

1 N 

(^n,y n ) n>no = ~y y^O^'lrn Vj^n) n>n ■ (5) 
i=l 

is another useful measure of the collective dynamics of the entire net- 
work, in particular in the context of the synchronization and stability 
of the CMS. 

• Time-averaged emergent orbit (t.a.e.o.), defined for each node as 

1 n 

(x[i],y[{\) = lim z2{x[i]k,y[i)k), (6) 

twoo n — uq * — ' 

k=no 

reduces the whole emergent orbit of a node to a single point in the 
phase space, that qualitatively captures the emergent motion. 

In the following we always refer to a particular type of orbit used. 

Structure of Emergent Orbits of Coupled Maps. As stated above, we 
consider CMS defined by Eq. ([2]) with fixed e = 0.9 at each node and the 
time-delay between interacting nodes r = 1, and vary the coupling strength 
between nodes Here we concentrate on the dynamical orbits of a selected 
node on the 4-star motif and a node on the tree graph. 

In Fig. [2] we show examples of three typical e.o. obtained for different 
fi -values for a node on the the 4-star motif. In particular, at weak coupling 
between nodes, u = 0.005, we have a chaotic motion with orbits wandering 
in the phase space in an irregular fashion (Fig.[2h&;d), resembling the orbits 
of the uncoupled standard map Eq. (TTJ. However, above a threshold coupling 
fi ~ 0.02 the CMS exhibits a regular motion. Orbits of different periodicity 
can be found in this region, such as the one in Fig. [2b. As a rule, the 2- 
periodic orbits are the ones which survive at large couplings u > 0.08. In 
the region between 0.02 < /i < 0.08, depending on the initial conditions, we 
also find a non-periodic motion with the fractal attractors, as the example 
shown in Fig[2b. The time dependence of the momentum coordinate in the 
respective three dynamical regimes is also illustrated in Fig.[2H-f. Clearly, 
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Figure 2: Three typical emergent orbits (a-c) and the respective time- 
dependence of the ^/-coordinate (d-f) for a node on the 4-star motif: (a,d) 
chaotic orbit at \i = 0.005, and (b,e) periodic orbit at [i — 0.02, both for the 
central node, and (c,f) strange attractor of a branch node at = 0.049. 



our CMS of Eq. (j2j) exhibits a variety of emergent dynamical orbits, in which 
both sensitivity to mutual couplings and to the initial conditions play a 
role. The appearance of the strange attractors with non-chaotic behaviors 
is a remarkable feature of the collective dynamics of our CMS, to which we 
devote the Section I4.2I Here we describe the differences between these three 
dynamical regimes in a quantitative manner. 

Dynamical Localization. In order to understand the nature of our CMS, 
we first analyze the chaotic e.o. occurring for arbitrarily small but non-zero 
couplings. Although these are chaotic orbits (see also the stability analysis in 
the Section I4~2"|) . in the presence of coupling they tend to localize in the action 
y— coordinate, as shown in Fig.[2h&d. Specifically, as opposed to the chaotic 
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diffusion in y-coordinate, which is well-known for the uncoupled standard 
map Eq. (pQ) (see e.g. [25]). the interaction between nodes in Eq. ([2]) inhibits 
the diffusion. Quantitatively, the mean-square distance < y 2 {t) > as function 
of time t for different coupling strengths \x is shown in Fig.[3^i. Clearly, the 
normal diffusion (7 ~ 1) observed for uncoupled maps fi = is replaced 
by the asymptotically localized orbits (finite m.s.r.) when [i > is turned 
on. A similar effect was demonstrated in the spreading of the band widths 
occupied by the e.o., Fig.[3b. The band width in y— coordinate increases as 
a power of time for uncoupled standard maps, whereas the stretching rate 
decreases with a non-zero coupling. Apart from quantitative differences, such 
localization occurs both on large trees and on small motifs. 




Figure 3: Dynamical localization of the emergent orbits for various coupling 
strengths: y-coordinate mean square distance (a) and average stretching of 
the band width (b). Plots obtained by averaging over an ensamble of e.o. of 
4-star branch node. Transients were also included. 



The observed dynamical localization is a clear argument that the collec- 
tive effects due to node-interaction are occurring. As we shall see in the 
following paragraphs, further decrease in the stretching rate of the band 
widths occupied by an e.o. triggers the appearance of a regular motion. 

Generalized Synchronization and Periodic Orbits. The appearance of pe- 
riodic e.o. is one of the central collective feature exhibited by the CMS 
Eq. ([2]). The trajectories, similar to the one in Fig.[2b&e, oscillates in the 
phase space between two localized groups of points in a way to maintain a 
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constant average y. For further quantitative study we use the time-averaged 
e.o. (x, y) [i] defined in Eq. (EJ) and consider the star's branch-node. As t.a.e.o. 
is represented by a point in 2D node's phase space, here we plot only its In- 
coordinate as function of the coupling strength \i , as shown in Fig.H^i. The 
color code represents (on log-scale) number of initial conditions that end up 
with a given value of y. 




Figure 4: (a) For a star-branch node: 2D distribution histogram of time- 
averaged emergent orbits y against coupling strength \x. Color-code: log- 
arithm of the fraction of e.o. corresponding to a given y. (b) Fraction of 
non-periodic e.o., fr(np), as function of fi for a star's branch node and an 
outer node of the tree. Periods up to 10000 are considered. 

The Gaussian distribution of y-values for small /x reduces to an organized 
set of clusters as the coupling strength is increased above \x ~ 0.01. With 
the disappearance of the chaotic e.o., classes of periodic e.o. appear. In 
this region, depending on the initial conditions, an orbit ends up with a 
given y- value, i.e., at one of the horizontal lines in FigHh. The e.o. perform 
an oscillatory motion around a fixed set of phase space points, which are 
located symmetrically around the respective average value y. The emergent 
periodicity is studied in more details later. A similar structure of orbits is 
found for each node on the graph, with robust patterns of groups of nodes 
having the same y value. Spatial patterns of such clusters of nodes depend on 
the local structure of the graph [20]. On the scale-free tree graph the number 
of possible groups gradually decreases with increased coupling strength and 
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eventually reduces to a single group centered around y — for fi ~ 0.07. 

Appearance of Non-Periodic Orbits. Furthermore, in the same region of 
couplings we find a scalar fraction of initial conditions leading to non-periodic 
orbits. In Fig.Hb we show the fraction of non-periodic e.o. against coupling 
strength both for a node on the 4-star and for a node on the scale-free tree. 
Nature of the non-periodic orbits in our CMS will be studied in more details 
in the next Section. Here we point out a systematic correlation between two 
curves in Fig.Hb, suggesting the importance of the dynamics at small-scale, 
4-star motif, in the genesis of the non-periodicity at the tree graph. 

3.1 Statistical Features of Emergent Orbits 

As the above results and Fig.H] suggest, the coupled chaotic maps defined 
by Eqs. (JTH2J) on tree graphs exhibit several types of collective dynamic be- 
haviors, apart from synchronized regular motion, which dominates at large 
couplings. In particular, we find (i) weakly localized chaos, {it) emergent pe- 
riodicity, and (m) non-periodic orbits with non-chaotic features (see later). 
With the appropriate statistical analysis we further differentiate between 
these dynamical regimes. 

Return-times statistics. The first-passage (or return-time) statistics rep- 
resents an interesting statistical measure of correlations in the dynamical 
systems. In this context, we monitor return of the trajectories to a marked 
small area in the phase space. The results, integrated over the phase space 
subset (x, y) G [0, 1] x [—50, 50] and averaged over many initial conditions for 
three representative values of the coupling \x are shown in Fig.(5k. Clearly, 
the distributions of return-times have different tails for the chaotic motion 
(the case with /i = 0.005), periodic (for fi = 0.021) and non-periodic (the case 
with [i = 0.048) behavior. In particular, the cooperative behavior leading to 
the power-law tail of the distribution in the case for \i = 0.048, suggests that 
the observed non-periodic orbits in this regime are non-chaotic. It is interest- 
ing that the differences between these dynamical behaviors are also expressed 
in the transient regime (approximately first 10 3 iterations in Fig.Hb). 

Distribution of Periods. Another measure of the cooperative behavior in 
the CMS is the emergent periodicity of the orbits at the level of an individual 
node coupled to the system and in the system-averaged trajectory. In Fig.Hb 
we show the distribution of periods tt of the periodic orbits for the tree- 
average trajectory and, in the inset, ranking statistics of the periods observed 
at individual nodes coupled on a tree and on the 4-star. It is interesting 
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Figure 5: (a) Distribution of return-times to a selected part of the phase 
space in three dynamical regimes corresponding to three different coupling 
strengths fi for star's branch-node. Phase space grid of 1000 x 100000 cov- 
ering (x, y) G [0, 1] x [—50, 50] was used and distributions averaged over an 
ensamble of initial conditions, (b) Distribution of periods of the network- 
averaged orbits on the tree for \x = 0.08. Inset: Ranking statistics of periods 
for a branch node on the 4-star and an outer node on the tree for fi = 0.21. 



that at the tree-like topologies the odd and even periodicities appear to have 
different statistics, which is especially pronounced at lower values of 7r and at 
weaker couplings. Although the large periods can be found at nodes on small 
graphs, our findings suggest that the large-scale structures are more efficient 
in building-up the power-law dependencies in the distribution of periodicity. 
The network-averaged trajectories, defined by Eq. (jSJ), exhibit the emergent 
periodicities with the power-law distribution shown in Fig.Hb. The curve 
can be well fitted by the q— exponential form [27] 

P{X) = B q [l-{l-q)^Y\ (7) 

with the characteristic period X = 63 above which the power-law tail ap- 
pears with the slope 1.66. The fit is compatible with the value of the param- 
eter q = 1.6. Often large deviations of the parameter q from unity can be 
related with a non-ergodicity of the underlying phase space. 
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4 Dynamic Stability and Evidence of SNA 



The analysis in previous Sections revealed that our CMS maintains differ- 
ent types of cooperative dynamical behaviors, depending on the coupling 
strength. In order to further substantiate the nature of the emergent dynam- 
ics, in particular the occurrence of the non-periodic orbits, here we discuss 
the appropriate dynamic stability of the system. 

We analyze the trajectory divergence [28] within small neighborhoods 
of the emergent orbit. As it was shown above, for the dissipative CMS of 
Eq. ([2]) the trajectories at long times tend to localize at a finite distance in 
the action coordinate. This localization implies that only limited divergence 
between different trajectories can be observed at t — > oo. Therefore, the 
appropriate stability measure is the initial divergence rate, which lasts for a 
finite time interval and determines the finite-time Lyapunov exponent [28] . 

4.1 Finite-Time Lyapunov Stability 

The Maximal Lyapunov Exponents (MLE) for a phase space point is gener- 
ally defined as the maximum rate of the divergence among all the trajectories 
starting in the vicinity of the considered point. In Fig. [6^ we show several 
examples of different initial divergence rates for different trajectories. As ex- 
pected, all curves reach a saturation (slope zero) at long times. However, the 
trajectory separation increases over certain number of initial iterations n with 
a well defined slope, before the saturation occurs. These slopes determine 
the finite-time MLE, denoted as A^ aa .. 

More precisely, the finite-time Lyapunov exponent A^ ax (x ) associated 
with the point x = (x , yo)[i] on a trajectory is defined by 



where M stays for a small neighborhood around the point x , d(., .) denotes 
distance and U n stands for the discrete time-evolution dynamics given by 
CMS Eqs. flUO]). Note that n should be appropriately selected for each tra- 
jectory separately. As before, we consider trajectories of a selected node 
coupled on the 4-star and on the tree. 

The procedure employed for computing the A^ a3 ,(x ) is the following: 

1. take an initial condition x (for all the nodes) and compute their e.o.; 
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Figure 6: (a) Behavior of the orbit divergence rates for four typical e.o. 
with different couplings /i . (b) 2D distribution histogram of A^ a:r for star's 
branch-node versus [i . Color-code reports the logarithm of the fraction of 
orbits leading to the same X max for a given n -value. 



2. focus on the final state (point) of a selected node and consider a ran- 
dom point in its close neighborhood, with the distance d defined as 
"Manhattan distance" d((x, y), (x , y )) = \x — x \ + \y — y \; 

3. iterate the dynamics for both points with systematically recording the 
rate that their distance ratio Eq. ([8]) changes in time and determine 
the optimal n; 

4. approximate the slope over the optimal estimation time n. Repeat the 
same procedure for a few other points in the same neighborhood and 
determine the maximal slope obtained this way. 

Clearly, A^ aa ,(xo) corresponds to only one point on the e.o., implying that, 
using different points, a spectrum of A^ ax -s can be determined to characterize 
the entire orbit. Furthermore, we define an average value \ max over the 
distribution of different Aj nax (x ) for a given orbit as: 

Knax — < ^Laxi^o) >x ee.O. (9) 

We compute the average \ max for a branch-node on the 4-star for many 
trajectories and different coupling strengths The results are shown in 
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Fig. [2b in the form of 2D histogram for /j, in the interval [0.0,0.08]. Not 
surprisingly, fairly large positive values of A^ aa , at fi < 0.01 are in a full 
accordance with the chaotic motion of weakly coupled nodes. More interest- 
ing is the appearance of the orbits with A^ a;c > in the range of couplings 
0.01 < ji < 0.08, where most of the trajectories are localized (cf. Fig.H^&b). 
Note, however, that the measured values of A^^. in this region are far smaller 
compared to the chaotic e.o. at small fi . This indicates that the nature of 
these e.o. might be fundamentally different. It should be also noted that, in 
this region of couplings we can find both positive and negative A^ aa , for differ- 
ent initial conditions (or different parts of the trajectory), as discussed later. 
Qualitatively similar picture is found for the trajectories of an outer-node on 
the large tree. 




0.08 




0.04 0.06 0.08 



(a) (b) 

Figure 7: Comparison of the fraction of orbits with a positive A^ ax and 
fraction of non-periodic e.o. for (a) a branch-node on the 4-star motif and 
(b) an outer node on the large scale-free tree. 



Next we compare the occurrence of the trajectories with A^ a;E > and 
the non-periodic orbits, mentioned before. Note that the periodicity is usu- 
ally determined after long times when the trajectory is localized, whereas the 
value of A^ ax is the average computed along the entire trajectory. Neverthe- 
less, as shown in Fig.[7h&b, we find a significant overlap in the fraction of 
the non-periodic orbits and the orbits with a positive finite-time Lyapunov 
exponent A^ ax . In the case of a 4-star's node the overlap is nearly complete, 
whereas in the case of a tree's node the qualitative similarities are clear. 
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4.2 Appearance of Strange Attractors 



The observed excessive number of non-periodic orbits, as compared to the 
orbits with an average A^ aa , > 0, suggests that we take a close-up look 
at the properties along the non-periodic orbits. The two examples of such 
orbits found for the branch-node on the 4-star for [i = 0.048 end-up at the 
strange attractors, parts of which are shown in Fig.[H^&b. A systematic study 
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Figure 8: Examples of strange attractors with (a) A^ ax = 0.118 and (b) 
Knax = —0.005, both for a branch-node on the 4-star. (c) and (d) Spectrum 
of A^ a;r (xo) values along the trajectories shown on (a) and (b), respectively. 



of the exponent A^ ax (x ), defined in Eq. ([8]), at different points x along 
the trajectory gives the spectrum of A^ a:r values. For the two trajectories 
in Fig.[8^&b the respective spectra are given in Fig.[8b&d. Note that the 
A^ ax , reported above, roughly corresponds to the peak value of the spectrum. 
The fractal attractors, as those shown in Fig.[8k&b, have a common feature: 
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the spectra have tails on both positive and negative values of A^ a3 ,(x ). 
Specifically, the attractor in Fig.[Sb has a negative distribution average value 
A^ ax close to zero. Attractors of this type are referred to as strange non- 
chaotic attractors SNA [28], indicating the atypical presence of a negative 
Lyapunov exponent accompanied with the fractal structure. 

We find a large number of strange attractors in our CMS. An incomplete 
search for fractal attractors gives the distribution of the peak values )\ max 
of the respective spectra shown in Fig.(9k, indicating roughly four different 
groups of such attractors in our CMS. 




max 




(b) 



Figure 9: (a) Distribution of }\ max for 10 3 trajectories of the branch-node on 
the 4-star for \i = 0.048. Averaging was done over spectrum with 100 values 
of A^ aa ,(xo) computed along each trajectory, (b) Distribution of return times 
for two trajectories, parts of which are shown in Fig.[8k&b, with smaller and 
large slope, respectively. 



Further interesting question is how the appearance of such attractors 
affects the collective motion of the CMS. Here we present the distributions 
of the return times to a given portion of the phase space, P(At), defined 
above. For the two attractors in Fig.[8ti&b the results for these distributions 
are shown in Fig. [9b. In both cases the distributions of return times exhibit a 
power-law tail, which can be fitted with the q-exponential function in Eq. ([7]), 
suggesting a non-ergodic dynamics near these attractors. Although both 
attractors in Fig.[8k&;b have quantitatively similar fractal dimension, close 
to 1.5, obviously their other features affect the return-time statistics. In 
particular, the larger slope and smaller characteristic scale X = 80 was 
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found in the case of the SNA of Fig.[Sb, whereas, at the attractor Fig.[Hk 
with a positive A^ aa ,, we find X = 8000 and an increased probability of long 
return times. Consequently, the parameter q in the expression Eq. (J7|), which 
is a measure of non-ergodicity, appears to be different, i.e., q = 1.6 for the 
attractor in Fig.[HK and q = 1.5 for SNA in Fig.[8b. 

5 Conclusion 

We have demonstrated that a variety of new dynamical effects occur in two- 
dimensional standard map with a non-symplectic coupling and time delay 
between coupled units on tree-like graphs. By focusing on a single node 
coupled to the tree, we were able to apply the established methods of the 
discrete-time dynamics, and thus closely follow how a coupled unit adapts to 
the evolution of the extended dynamical system. In particular, as a conse- 
quence of the coupling we find the dynamical localization of orbits, emergent 
periodicity, non-periodic orbits and occurrence of fractal attractors in a wide 
range of coupling strengths. The self-organized dynamical behavior of mutu- 
ally coupled maps in our CMS is characterized by power-law dependencies, 
specifically in the emergent (tree averaged) periodicity and in the return time 
intervals near the strange attractors. 

We have shown that some of the strange attractors in our CMS have non- 
positive Lyapunov exponents, similarly to the strange non-chaotic attractors 
[28l [29l [30] found in the quasi-periodically forced maps. It should be stressed 
that in our two-dimensional coupled maps no external driving was imposed. 
Contrary to the forced isolated maps, in our CMS different coupled units on 
the 4-star or on the scale-free tree, provide the dynamical input to each other. 
Full understanding of the routes to the strange non-chaotic attractors in this 
coupled-map system requires additional study. We speculate that both the 
observed localization of orbits and time delay between coupled units, which 
pervents full synchronization of neighboring nodes, play an important role. 

With the parallel analysis of the dynamics at the 4-star motif and on the 
large tree, in this work we obtained certain quantitative arguments, which 
imply how the self-organized dynamic on the large scale networks occurs. 
Specifically, the 4-star appears as a basic dynamical unit on which the non- 
periodic orbits and strange attractors appear and spread on the entire tree, 
cf. Fig.Hl The role of cycles on the graph as well as the effects of dif- 
ferent/variable delays between units are left for a separate work. Further 
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interesting subjects in the context of our CMS, as the appearance of partial 
synchronization patterns (see also [20J) and the problem of coexisting attrac- 
tors on the graph, are left for the future study. 
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